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In this article we consider a piston modelled by a potential in the presence of extra dimen- 
sions. We analyze the functional determinant and the Casimir effect for this configuration. 
In order to compute the determinant and Casimir force we employ the zeta function scheme. 
Essentially, the computation reduces to the analysis of the zeta function associated with a 
scalar field living on an interval [0, L\ in a background potential. Although, as a model for a 



piston, it seems reasonable to assume a potential having compact support within [0, L], we 
provide a formalism that can be applied to any sufficiently smooth potential. 

I. INTRODUCTION 

In recent years piston configurations have received a surging interest in the Casimir effect 
community. The main reason for this fact is that pistons allow for an unambiguous prediction of 



forces which turn out to be divergence free 



__]. The piston is usually represented by an infinitely 



thin movable plate at which the field has to satisfy some ideal boundary conditions. Different 
boundary conditions and various shapes of cross-sections have been analyzed and, as expected, the 
force heavily depends on the different possible choices, see, e.g., , llal . It is the aim 

of this article to represent pistons of finite thickness by compactly supported potentials. Physical 
properties of the pistons are encoded in the spectrum of the ordinary differential operator 

P:= -^ + y(x) ' (L1) 

where x € I = [0, L] and V(x) models the piston. The points x = and x = L represent the 
positions of the fixed plates with suitable boundary conditions chosen. The piston at position a is 
represented by the potential V(x) which has a support strongly concentrated around a. 

It is this problem to which an analysis in the space I x J\f reduces after separation of variables, 
where N describes the cross section of the piston and additional dimensions, represented by a 
smooth Riemannian manifold possibly with boundary. 
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Although ultimately, in the context of pistons, our interest is in potentials with compact support 
within /, the formalism is developed for general potential V(x). We will use the zeta function 
scheme to evaluate various quantities of interest; for introductions to spectral zeta functions and its 
app.ica.ioas in plly sics see JlQQ. The relevant ,e t a fusion is _ t ed as a c„n t „n r i^a! 
where the boundary value of the unique solution to a given initial value problem enters. This 
representation is briefly described in Section |TJ and subsequently used to evaluate the functional 
determinant, the Casimir energy and force. In order to analyze these quantities, the analytic 
continuation of the zeta function needs to be constructed. This will entail the knowledge of a 
certain asymptotic behavior of solutions to initial value problems, which is obtained in Section 
IIIII using standard WKB techniques can be used to find the 

functional determinant, Section IIV| and the Casimir energy and force, Section [V] In Section IVII 
we restrict to compactly supported potentials and describe the resulting simplifications. In the 
Conclusions we summarize the most important findings and outline possible further applications 
of our approach. 

II. CONTOUR REPRESENTATION OF THE PISTON ZETA FUNCTION 

Let M = [0, L] x J\f, where J\f represents the cross section of the piston and the additional 
Kaluza-Klein dimensions. For simplicity we assume Dirichlet plates at x = and x = L and we 
assume a sufficiently smooth potential V{x) depending only on x G [0, L\. With y £ J\f, the relevant 
energy eigenvalues for a scalar field are determined by the second order differential operator 



together with Dirichlet boundary conditions at x = and x = L and unspecified boundary condi- 
tions at the boundary of N . Using separation of variables we write eigenfunctions in the form 





4>{x,y) = X(x)tp(y) 



where the <p(y) are assumed as eigenfunctions of the Laplacian on A/", 




This implies that the eigenvalues A of L are given as 

A = v 2 + rjl 
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where v 2 is determined by 



d 2 \ 

-^2 + V(x)j X„(x) = v 2 X v {x), X u (0) = X U (L) = 0. (2.3) 
Formally, the zeta function of L can therefore be written as 

C(s) = Yy + V 2 e r s for^ S >^, (2.4) 

with D = dim(M), the dimension of M. Note that without specifying J\f the spectrum r\n is not 
known, and also v 2 will not be known unless V(x) is one of the very few potentials allowing for a 
closed solution of eq. (|2.3p . 

Despite this lack of knowledge an analytical continuation of £(s) in eq. (j2.4j) can be constructed 
and properties of on M can be given in terms of the zeta function of N defined by 

Ca^) = E^ 2S forK S >^i. 

I 

The strategy for the analysis of in eq. (|2.4p is to rewrite the series as a contour integral using 

. Instead of considering the eigenvalue 



10, 



m 



the argument principle or Cauchy's residue theorem 
problem in eq. (|2.3p we consider the initial value problem 

d 2 \ 

^-2 + V[x)\ u^x) = fi 2 u^x), u„{0) = 0, u^(0) = 1, 

where fj, Ed. The eigenvalues v 2 of the original problem are recovered as solutions to the secular 
equation 

M L ) = °; ( 2 - 5 ) 

note, that u^(L) is an analytic function of [i. Eq. (|2.4p . for !Rs > D/2, can then be rewritten as 

^ S ) = ^E / d^ 2 +r, 2 r S j-^u,{L), 
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where 7 encloses all solutions to eq. (|2.5p . which are assumed to be on the positive real axis; the 
changes necessary when zero modes or finitely many negative eigenvalues are present are given in 



Let us next consider the contributions from each I by analyzing 

1 I" - ' .2 , 2n- s d 



0(s) = ^— J dfi(n z + rif) s — lnn M (L). 
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Deforming the contour, as usual, to the imaginary axis we find 



oo 

(l(s)= S -^ [ dk(k 2 - V jr s ^-lnu ik (L), (2.6) 
7r J ak 

valid for 1/2 < < 1. In order to construct a representation of Q(s) that is valid in a region 
3ts < 1/2, as is needed for the functional determinant and the Casimir energy, we add and subtract 
the large- k asymptotics of Ui k (L). 

III. ASYMPTOTIC BEHAVIOR OF BOUNDARY VALUES FOR DIFFERENTIAL 

EQUATIONS 

The next mathematical task therefore is to determine the large-/c asymptotics of the unique 
solution to the initial value problem 

d 2 \ 

^ + V(x) + k 2 ju ik {x) = 0, u ik (0) = 0, u' ifc (0) = l. (3.1) 

We note that the differential equation in (|3.ip has exponentially growing and exponentially decaying 
terms [l[ 3]. Although, ultimately, for (j2.6|) we will only need the exponentially growing part, 
we first have to consider linear combinations of the two in order to be able to impose the initial 
conditions in eq. (|3.ip . 

It is convenient and standard to introduce 

S(x,k) = d x lmp k (x), 

where ili k {x) satisfies 

.— + V{x) + k 2 U k {x) = 0. (3.2) 

The differential equation satisfied by S(x, k) turns out to be 

S'(x,k) = k 2 + V(x) - S 2 (x,k), (3.3) 

where the prime indicates differentiation with respect to x. As k — > oo, the function S(x,k) can 
be seen to have the asymptotic form 



i=-l 
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where the asymptotic orders Si(x) are given by 



= ±1, S (x) = 0, S 1 (x) = 



(3.4) 



It is clear that an arbitrary number of asymptotic orders can be easily evaluated using an algebraic 
computer program. The two different signs in (|3.4p produce the indicated exponentially growing 
and decaying solutions ipk( x ) °f (|3.2[) . We denote solutions of fj3.3[) corresponding to these two 
signs by S + (x, k) and S~(x, k). The associated solutions of (|3.2p then have the form 

^±( x ) = ^±exp | J dt S ± (t,k) 
The large-A; behavior for un~(x) is obtained by considering the linear combination 

u ik (x) = ,4+ exp < 



' X 




' x 1 


< J dt S + (t,k) 


> + A exp < 


J dt S~(t,k) 


.0 




,0 J 



together with the initial conditions in (|3.ip still to be imposed. These initial conditions imply 

1 



= —A' 



.4 H 



S+(0,A!)-5-(0,jfe)' 

note, that without including the S~(x, k) part the initial conditions could not be satisfied. 

We are now in the position to write out the large-fc behavior for Ujfc(L). Let E(k) denote 
exponentially damped terms as k — > oo. First, we see that, as k — > oo, 

( L 



Uik(L) 



1 



S+(0,k) - S~(0,k) 



exp < 



J dt S + (t,k) I +E(k), 



K 



and therefore 



lnu ik (L) = -\n(S + (0,k) -S~(0,k)) + J dt S + (t,k) + E(k) 

o 

oo 

= -\n{2k) + kL + s ^ j d j k- j + E(k), 

j=0 

where the dj's are easily determined from eq. (13.4p . For instance, the first six are given explicitly 
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by 



d 



0, d^lJdtV®, d 2 = ~[V(L) + V(0)], 
o 

L 

ds = \[V'(L)-V'm- l -J dtV 2 (t), 



C?4 



~[V"(L) + V "m+ l -[V\L) -V\0)], 



db = h [v{ ' ){L) " y(3)(0)] " U V{L)V ' {L) " v ^ v '^ + rj dt v ' {t) 



L 

~ h I dt v ^ v "^- 



(3.5) 



In what follows, the potential is always assumed to be as smooth as necessary for the asymptotic 
orders given by these formulas, and higher ones if needed, to be well defined. 

Subtracting and adding the asymptotic behavior up to the order k~ N , the zeta function naturally 
splits into two parts, 

where 



oo 

<**(.) =^/* 

oo 

c { r\s) = s ^ fdk (k 



N 



lnUjfc(L) - kL + ln(2k) - ^djk' 



3=0 



N 



: nf ) s 



Ve 



kL-ln{2k) + J2djk- j J- . 

3=0 



(3.6) 
(3.7) 



The /c-integrals in (i as \s) are easily done, yielding 



* r s + 



After summing over £, the representation obtained is then valid for 1 > Sfts > (D — N — 2)/2 and 



7 



it reads 

00 ( N 

C (/) (*) = ^Ej dk (**" ^' S ^ - *L + In(2fc) - £ § } , (3.8) 



'/( 



j=0 



N y(s + 2-\ 



j= i i^tj 

In particular, choosing N = D — 1, respectively N = D, the representation can be used to compute 
the determinant, respectively the Casimir energy. This will be done in the next sections. 

IV. FUNCTIONAL DETERMINANTS 

In this section we will evaluate the functional determinant, or, equivalently, C'(0), using the 
representation of £(s) given by eqs. (|3.8p and (|3.9p . The contribution from C ( s ) i s trivially 
obtained and it reads 



c^'(o) = -E 



lnu i% (L) - Lrj£ + ln(2%) - ^ dj7] e 

3=1 



as the sum over £ converges by construction. For the evaluation of the contribution from C^ as \s), 
let us note that in the general situation considered, namely J\f can be a manifold of any dimension 
with or without boundary, for j = 1, (D — 1), we have the following Laurent expansion 

CM f ~ + A = ^Res Cv(j/2) + FP Cv0'/2) + O(e). (4-1) 

Using standard properties of T-functions, this immediately gives 

C<">'(0) = -L ■ FP Cjj + L ■ Res Or (-i) (-2 + In 4) - i&(0) 

- g d, (FP & (0 + Res <y (0 ( 7 + i, (0 ) ) . 

Adding up these two pieces gives the answer for the functional determinant on M in terms of the 
spectral zeta function on N '. Once the 'base' manifold M is specified, more explicit results can be 
given. 
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V. CASIMIR ENERGY 



For the Casimir energy we set N = D and use again the Laurent series (|4.ip for Q/(j/2 + e). 
In this case, we find 

D 



C (/)(_l /2 ) = fdki^-^^llnu^-kL + ^k)-^ 



hi 



C M(_l/ 2 + e ) = - i ^Cat(-I) - -Res Caa(-1/2) 



1 I ^ 



e 47r 



^iCAr(O) 
2vr 



+^ (CV(-l) + Cvr(-i) (-1 + in 4)) - Ifp Cat(-i/2) 

+^(CV(0) + C^(0)(-2 + ln4)) 



D 4 r ft 1 



,=2 2 ^ r(| 



with the harmonic numbers 



fpCat 



R^kv(^) [%-3)/2-2 + ln ! 1 ■ 



- 1 



i=i 



Multiplying by 1/2 and adding up the two pieces, the Casimir energy follows. Again, when M is 
specified more explicit results can be given. In general, the Casimir energy is divergent, but as 
seen below, the resulting forces on pistons are finite. 



VI. COMPACTLY SUPPORTED POTENTIALS 



In order to reasonably talk about the force on a piston modeled by a potential we now assume 
the potential V(x) to have compact support within the interval [0, L\\ namely, we assume that 
it does not vanish for x S [a — e, a + e] C [0, L\. This can be considered a model for a piston of 
thickness 2e. In this case the asymptotic behavior of Uik(L) will be independent of a as the integrals 
over V(x) and its powers and derivatives are independent of a (as long as the support is within the 
interval [0, L]). This can be seen explicitly in (13.5|) . It is here that sufficiently smooth potentials 
become a necessary assumption in order for these formulas, and the corresponding ones for higher 
orders, to be well defined. These formulas also simplify further because V(0) = V(L) = 0. The 
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corresponding result can be used to write down the Casimir energy. Because of the independence 
of the asymptotic behavior of Uik(L) from a, in these circumstances we immediately obtain 

oo 

The force, in particular its sign, is encoded in boundary values of an initial value problem to an 
ordinary differential equation. 

Results for the determinant are also easily written down from Section IIVI with the simpler dj 's 
used. 

VII. CONCLUSIONS AND OUTLOOK 

In this article we have provided a formalism that allows for the evaluation of functional de- 
terminants and Casimir energies and forces for the configuration of a generalized piston. Results 
are as explicit as they can be without specifying the cross section and the additional Kaluza-Klein 
dimensions that might be present. 

With these results available, one can obtain very explicit answers for a given cross section of 
the piston and specified Kaluza-Klein dimensions. Furthermore one can use the potential to mimic 
material properties of the piston. For all cases, with very few exceptions for particular potentials, 
a numerical determination of the boundary value Ujfc(-L) will be necessary. 

Along the same lines different boundary conditions at x = and x = L can be considered. 

Work along these lines is currently in progress. 
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